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Abstract 

We use the dynamical mean field method to study a model of electrons 
Jahn- Teller coupled to localized classical oscillators and ferromagneti- 
cally coupled to "core spins", which, we argue, contains the essential 
physics of the "colossal magnetoresistance" manganites Rei_xAxMn03. 
We determine the different regimes of the model and present results 
for the temperature and frequency dependence of the conductivity, the 
electron spectral function and the root mean square lattice parameter 
fluctuations. We compare our results to data, and give a qualitative dis- 
cussion of important physics not included in the calculation. Extensive 
use is made of results from a companion paper titled: "On the Fermi 
Liquid to Polaron Crossover I: General Results" . 
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I. INTRODUCTION 



The doped rare earth manganites have been studied for many years |l|,0] and 
interest in the materials has revived following the recent discovery of extremely 
large magnetoresistance in some members of the family [^. The chemical formula 
is Rci-xAxMnOa, with Re a rare earth element such as La or Nd and A a divalent 
metal ion such as Sr or Ca. The electronically active orbitals are believed to be the Mn 
d-orbitals ||T],3 and the mean d-occupancy is 4— x. Each Mn ion feels an approximately 
cubic crystal field which splits the Mn d-levels into a t2g triplet and an Cg doublet 
0]. The t2g levels are believed to lie substantially (~ 5eV) below the Cg levels. 
On-site Coulomb interactions are apparently strong enough that no d-orbital may be 
occupied by more than one electron. Further, all electron spins in Mn-d orbitals are 
aligned by a large ferromagnetic Hunds-rule coupling. The Coulomb and Hunds-rule 
interaction energies have not been measured directly, but there is substantial indirect 
evidence that they are large. For example, at 0.2 < x < 0.5 (precise values depend 
on Re and A) the ground state is ferromagnetic, and the observed magnetization is 
consistent with all 4 — x electrons on each Mn being lined up in the maximal spin 
state 1^, suggesting a large Hunds coupling. Also, ReMnOa undergoes a structural 
phase transition at T 800K, which has been shown |^,|^ to be due to a staggered 
(tt, tt, tt) ordering of Jahn- Teller distortion of locally Cg symmetry. This would not 
occur unless the Cg orbital were singly occupied, which in turn implies that the t2g 
orbitals are also singly occupied, suggesting a large on-site Coulomb interaction. 

The resulting physical picture is that 3 of the (4 — x) d electrons fill up the t2g 
levels, forming an electrically inert core spin Sc of magnitude Sc = 3/2. The remaining 
1 — X electron goes into a linear combination of the Cg levels, and is coupled to Sc 
by a Hunds rule coupling Jh which is presumably large, but has not been directly 
measured. Okimoto et al. have recently presented an interpretation of optical data 
implying that JhSc ~ 1.2eV We shall argue below that their interpretation is not 
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correct and that JrSc is rather larger. Certainly, the conventional [|l| wisdom is that 
the limit JhSc ^ oo is appropriate, so one only need consider configurations with Cg 
electrons parallel to core spins. 

The Rei-xAxMnOs materials display a wide range of interesting physics. For 
< X < 0.2 (all x-values are approximate, and depend on Re and A) the materials 
are insulating at all temperatures and are antiferromagnetic or ferrimagnetic at low 
T. For 0.2 < X < 0.5 the low-T phase is a fully polarized ferromagnetic metal. As 
the temperature is increased for 0.2 < x < 0.5, there is a ferromagnet to paramagnet 
transition, which may be of first or second order, at a Tc(x) ~ 300K. In the param- 
agnetic state the material may be either "metallic" (in the sense that dp/dT > and 
P ^ PMott) or "insulating" (in the sense that dp/dT < and p > pMott)- (Here pMott, 
the Mott "maximum metallic resistivity", is about ~ 1000pi7-cm and corresponds 
to a mean free path of order pp^ 0). Insulating behavior occurs at lower x and 
metallic behavior at higher x. There is for all x a very pronounced drop in p as T is 
lowered through Tc, and in this regime the resistivity has a very strong magnetic field 
dependence. The "colossal" magnetoresistance of interest here occurs for x such that 
the material is in the "insulating" regime at T > T^ but is a metallic ferromagnet at 
T < Tc. Finally, at x > 0.5 the low T state is charge-ordered, antiferromagnetic and 
insulating |]I|,p!0|. We do not address the physics of this regime here. A qualitative 
phase diagram is shown in Fig. 1. 

Most models [p!T| -[T5| of electron transport in Rei_xAxMn03 have emphasized the 
"double exchange" phenomena caused by the large Hunds coupling Jh. The essence 
of double exchange is that when an electron hops from site i to site j it must also 
go from having its spin parallel to SJ, to having its spin parallel to S^,; the hopping 
amplitude, ty, thus depends upon relative spin orientation [|Tl|]. For two fixed sites i 
and j it is possible to choose phase factors so that ty — > + = ty cos [6*^/2]. 

The double exchange phenomenon gives an obvious connection between electron 
hopping and magnetic order: disorder in the spins implies randomness in ty, which 
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decreases below Tc or in a field. This effect seems very likely to be related to the 
"colossal" magnetoresistance observed near T^. However, two of us and P. B. Little- 
wood have recently argued that models involving only double exchange cannot explain 
the observed resistivity |I6[. The essential point is that in materials exhibiting "colos- 
sal" magnetoresistance the resistivity at T > T^ is much larger than the Mott limit 
and moreover rapidly increases as T decreases. Indeed, as shown in Appendix A, the 
observed T > Tc resistivities are so large that a classical description involving par- 
ticles incoherently hopping from site to site with a hopping probability W <^ keT/ft 
is appropriate. In models involving only double exchange the scattering produced 
by spin disorder is simply not large enough to cause such insulating behavior. A 
straightforward calculation 0,11^] shows that if the spins are completely decorre- 



lated one finds ppi ~ 3, i.e. W ~ ty ^ kT. More sophisticated arguments involving 
localization and phase factors are shown in Appendix B not to change this conclusion 
significantly. Therefore we believe some additional physics not included in the double 
exchange-only model must be important. This conclusion is not universally accepted 

One possible source of this extra physics is the "Hubbard U" effect of the on-site 
Coulomb interaction which produced the Hunds coupling in the first place. While this 
is undoubtedly quantitatively important, we do not believe it is the primary cause 
of the observed insulating behavior, essentially because away from commensurate 
densities (such as one electron per site) canonical Mott insulating materials such as 
the high T^. superconductors or other doped transition-metal oxides have resistivities 



which are rather less than the Mott limit and which decrease with temperature [|18| , p!9 
in stark contrast to the behavior observed at T > Tc in Rei-xAxMnOs. 

We proposed that the crucial additional physics is a strong electron-phonon 
coupling, which localizes the conduction electrons as polarons at T > Tc and smaller 
X, but is weakened in the T < Tc ferromagnetic state, restoring metallic behavior. 
We argued that this is possible because the behavior of electron-phonon model is 



controlled by a dimensionless coupling parameter which is the ratio of an interaction 
energy to the electron kinetic energy. The double exchange physics implies that ferro- 
magnetic order increases the electron kinetic energy, thereby decreasing the effective 
coupling strength. Also, a recent analysis [^] of the structural distortion observed 
1^ in LaMnOa showed that the electron-phonon coupling is indeed strong. 

In this paper we present a detailed study of a model of electrons coupled to core 
spins and to phonons which we believe confirms the importance of electron-phonon 
interactions. We use a "dynamical mean field" method which has previously been 
extensively applied to interacting problems without double exchange |2l| and has 



been used by Furukawa to study models involving only double exchange |]T5[. The 
model we study does not capture all of the physics of Rci-xAxMnOs; in particular. 
Coulomb effects and quantum and intersite terms in the phonon Hamiltonian are 
omitted and an oversimplified electron-phonon coupling is used. We therefore cannot 
quantitatively compare our results to experiment. The qualitative agreement we 
obtain is however compelling. 

Other workers have also studied electron-phonon effects in manganites and related 
materials. Emin, Hillery and Liu studied a theoretical model of a single bound po- 
laron coupled to spin waves and found a temperature dependence of the polaron size 



which they argued could be related to transport anomalies at Tc observed in EuO 
Their work, however, is based on models in which double exchange physics is not rele- 
vant, reflecting the different physics of EuO. Roder, Zang and Bishop used variational 
wavefunction techniques to examine the interplay between electron-phonon interac- 
tion and double exchange |2^. Their work is in a sense complementary to ours. They 



have incorporated quantum phonons and have presented some results on intersite 
phonon correlations, but their technique seems to work best at low temperatures. To 
calculate properties at and above Tc they resort to a dilute limit approximation which 
amounts to the study of a single carrier in a deformable medium. Also, they have not 
presented results for transport and optical quantities; their main result is a calcula- 



tion of the coupling dependence and doping dependence of Tc. Their results for the 
coupling dependence are very similar to ours; their results for the doping dependence 
are based on an assumption which, we argue, is not justified by the arguments in 
their paper A comparison of their results for Tc to ours is given in Section IV, and a 
further discussion is given in the conclusion. 

The rest of this paper is organized as follows. In Section II we define the model and 
the approximations used. In Section III we give the qualitative physics of the model, 
distinguishing different regimes and presenting the behavior of physical quantities in 
each. In Section IV we present the results of a detailed numerical study of the model 
at half-filling. In Section V we present and discuss results at other dopings. Section 
VI is a conclusion, in which the relation of results to data is analysed and the effects of 
omitted interactions are outlined. Appendix A discusses the theoretical interpretation 
of the observed resistivities, in particular the relation to the Mott minimum metallic 
conductivity and to transport by classical particles while Appendix B discusses in 
more detail the resistivity of the double exchange-only model. An announcement of 
some of the results of this paper has been submitted elsewhere . 



This paper relies heavily on results of a companion paper which uses the dynamical 
mean field method to study electron-phonon interactions in models without double 



exchange |25], to which we refer henceforth as I. 



II. MODEL AND APPROXIMATIONS 

We study a model of Mn Cg d-electrons coupled to core spins SJ, and phonons, 
with Hamiltonian 

H = Hband + Hd-ex + Hgl-ph + Hph. (1) 

Hband describes the hopping of d-electrons between sites i,j of a lattice. We take the 
electrons to have two-fold orbital degeneracy labelled by a Roman index (a,b) and 
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two fold spin degeneracy labelled by a Greek index Explicitly, 

Hband = - ty''djl,„djba - /iXl'^Jaodiaa- (2) 

(ij)abcj iaa 

The hopping matrix element t^*^ is a real symmetric matrix whose form depends on 
the choice of basis in ab space and the direction of the i — j bond. The precise form 
will not be important in what follows. 

The interaction of itinerant electrons with the core spins is given by 

Hd-ex = — Jh X! ■ dJaa(?a/3dia/3 (3) 
iaa/3 

As discussed below, we shall take the hmit JrSc oo. 

We assume a Jahn- Teller form for the electron-phonon coupling. Previous analysis 



||6|j20|| has shown that this coupling is strong in LaMnOs, so it may be expected to 
be strong also in doped compounds. Thus we focus on lattice distortions which split 
the on-site orbital degeneracy of the eg levels. Physically, these correspond to Cg 
symmetry distortions of the oxygen octahedra around an Mn site. Mathematically, 
we may parametrize a local eg distortion by a magnitude r and an angle 6, and define 
a two-component vector r = (rz,rx) with r^ = Tcos(f) and rx = rsin0. The coupling 
of this to the Cg levels is prescribed by group theory 0] to be 

Hel-ph = g E ■ ^L^'^diha. (4) 
iabcK 

Here f = (r^, r^) is a vector of Pauli matrices acting in orbital space. 

It is important to note that the coupling written in Eq. HI is not the only physically 



relevant one. In ref it was argued that a Mn site with no Cg electrons would 
induce a breathing distortion of the surrounding oxygen ions, and that this breathing 
distortion played an important role in determining the x dependence of the structural 
phase boundary. We have not included this coupling in the present calculations, but 
will qualitatively discuss its effects in the conclusion. 

In order to obtain a tractable model we assume Hph describes classical harmonic 
oscillators of stiffness k which are furthermore independent from site to site. Thus, 
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Hph = E^krf. (5) 

i 

Despite the simplifying approximations, the model defined by Eq. ^ is not solv- 
able except in certain limits. To obtain results, we adopt the "dynamical mean field" 
approximation, which becomes formally exact in a limit in which the spatial dimen- 
sionality d — > oo and has been shown to provide a reasonable description of models 
of interacting electrons in d = 3 Recently, the technique has also been applied 
to the double exchange-only model defined by Eq. |l] with Hd-ph = [l^ . 



The dynamical mean field method is based on an assumption about the electron 
Green function G^^(p,ci;). In general this is a tensor in spin and orbital space which 
may be written 

G$(p, uj)] = ^ - + muE^^(p, u) (6) 

Here ep*^ is the dispersion implied by Eq. ^ and S is the self energy due in the present 
problem to Hei_ph and H^-cx- The fundamental approximation of the dynamical 
mean field method is the neglect of the momentum (p) dependence of S. This is a 
reasonable approximation because models of the form of Eq. |I| (such as the usual 
Migdal-Eliashberg electron-phonon Hamiltonian) generally lead to a self energy with 
a weak momentum dependence in d = 3 ||26|. If the momentum dependence of S 
may be neglected, then all physical quantities may be expressed as functionals of the 
momentum-integrated Green function Gloc, given by 



G?oc««H= 7^,G%{p,u;) (7) 



'LOCa/3l^J - J (27r)3 

We shall assume that there is no long range order in orbital space, so Gloc and also 
S(ci;) must be proportional to the unit matrix in ab space. We shall allow for the 
possibility of ferromagnetic order, and shall take the ordered moment parallel to z. 
We may then write 

GfoCa;3M=goH+glM^^ (8) 



and a similar equation for S. We simplify Eq. |^ by first writing go,i = 
^TrabTrQ,/? / ^|^G^^((j^^)°'^, then introducing at each p the rotation matrix Rp*^ which 
diagonalizes e^*^, i.e. 



,ab 



si 
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Rp^ 



(9) 



and finally exploiting the cyclic invariance of the trace. We obtain 



(10) 



go(a;) = iTrabTra/3 / depV{ep)[uJ - Cp + /i - T.ap{uj)] ^ 
gi(u;) = iTi^^Tiafi J dep'D{ep)a^[uj - ep + fi - J:af3iuj)]~^ 

where ep is either of the eigenvalues of (they are related by symmetry operations) 
and V is the density of states, which we take to be semicircular with width D = 4t: 



V{ep) = JAt^ - e2/27rt^ 



(11) 



Because Gloc is momentum independent and involves two independent functions, it 
must be the Green function of some effective single site model involving two mean 
field functions ao and ai. This model is described by the partition function 



Zloc = / rdrd0 / d^ncexp[SLoc 



(12) 

Here r and (f) are the classical oscillator coordinates introduced above Eq. ^ 
= Sc/|Sc|, and the integrals are simple integrals rather than functional integrals 
because we have taken r,0 and Qc to be classical. 
The effective action Sloc is 



Sloc = -2 +Trln[ao(iu;n) + a-i{iujn)(^z + Jh^c ■ ^af3 



(13) 



Here we have added a term coupling the core spin to an external field hgxt- One 
could also couple the external field to the Cg electrons, but the factor of 3/(1 — x) in 
size of moment means that this coupling is unimportant. 
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The mean field parameters ao, ai in Eq. |T3] are determined as follows |2T[]. One 
obtains the local Green functions gQ^^ = |(51n[ZLoc]/<^ao,i, defines from these self 
energies So,i = ao_i — (gQ^^*^)"^, and demands that = Sq + SiO"z) reproduces gg?*^ 
when used in Eqs. |10] and ^. For the semicircular density of states the resulting 
equations may be written 

a„M =. + ^-$«|^fl, ^^^^ 

The factor of four is that appearing in Eq. |10|. 

These equations simplify in the "double exchange" limit JrSc oo. The argu- 
ment of the Tr In in Eq. |l^ is a matrix in the direct product of spin and orbital space. 
It has four eigenvalues, ao ± A ± gr with A = |aiz + JhSc^]. These are, of course, 
independent of the angle describing the phonon. For JhSc ^ t the eigenvalues at 
ao — A ± gr correspond to high energy states which do not affect low energy phenom- 
ena. Further, from Eq. it is clear that ai is of order t, so ai ^ JrSc and we may 
approximate A ^ JrSc + a-iQ^- We may then absorb the constant term JhSc into ao 
and fi, rescale the parameters and define the spin angle 6 via Qz = cos{6) obtaining 

x^ 

S(x, 0) = -7^ + T. ln[(bo + bicos^)' - Ax^] + hocos0/T. (15) 



Here x = ryk/t, bo,i = ao,i/t, A = g^/kt, ho = hcxtSc/t, and T, to and /i are 
measured in units of t. The mean field equations become 

bo =u + ^,-y- xdx/^ d(cosg)P(x, 6) 

bi = -i Xdx J\ d(cOS^)P(x, e)cOSe . (bo+bicosg) 



(16) 



(bo+bicose)2-Ax2 



with 



P(x,e) = -^exp[SLoc(x,^^)] (17) 

^LOC 

These equations differ from those discussed in I by the presence of the angular 
integral and by the quantity bi, which expresses the spin dependence of G. Expressions 
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for physical quantities are also slightly different from those used in I because we must 
keep track of the spin dependence. 

The momentum integrated Green functions has components parallel (t|) and 
antiparallel (|J.) to the magnetization. The off-diagonal (|J,) components vanish. We 
have 

G};i,c(^) = uJ + ^I- bo(^) - biH (18) 
G[|)c(^) = uj + fi- bo(^) + bi(^) (19) 

We shall be interested in the spectral function 

A(tu) = -TrImGLoc(t^ + ^5)/% (20) 

The conductivity is given by 

= |^y'depP(ep)T5]Tr[G(p,i^)G(p,ia; + in)] (21) 

where the factor of two comes from the trace over orbitals. Here G is a diagonal 
matrix in spin space and we have set e = t = 1. 

Another interesting quantity is the electron kinetic energy, K, defined by 

K = TrabTr,^ J -^ef {d^^^d.^p) . (22) 

By use of the relation between the expectation value and the electron Green 
function, of Eq. ^ with momentum independent self energy, and of the arguments 
leading from Eq. ^ to Eq. |10| and the mean field equations, we obtain 

K = 2T [G[l,c(^n)] ' + [Gi^ocM] ' . (23) 

n 

The magnetization m is given by 

m= / xdx / dcos^cos^P(x,^). (24) 
Jo J~i 

In these units the T = value of m = 1. 
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We shall also be interested in the mean-square lattice distortion, x^, given by 



We conclude this section by mentioning numerical methods. We use the proce- 
dures described in I, and handle the additional angular integral by a twenty-point 
Legendre formula. Computations are of course more time consuming because of the 
extra integral involved. We found it convenient first to locate the magnetic transition 
temperature Tc and then to perform calculations at T > Tc using equations obtained 
by forcing bi = 0. Convergence difficulties arise for temperatures near Tc; these are 
presumably related to critical slowing down near the magnetic phase transition. We 
found that an accurate value for Tc was most conveniently obtained by computing 
several values of m in the range 0.15 < m < 0.3 (0.02 < m^ < 0.1) and finding Tc by 
fitting to m^(T) = a(Tc — T) with a and Tc fit parameters. 

In previous work p6| , p^ we had also used an alternative method (which we termed 
the projection method) based on the observation that by choice of an appropriate local 
spin reference frame one may map the model into one of spinless fermions moving in 
a lattice with a spatially varying hopping determined by the local spin orientations. 
We further argued that within mean field theory one could approximate this hopping 



by t(m) = ^(1 + m2)/2, thereby simplifying the problem to one of spinless fermions, 
with hopping t(m), coupled to phonons. Finally we argued that one could construct 
a mean-field magnetic free energy by combining the m-dependence of the free energy 
of the auxiliary problem with the entropic term from the conventional mean field 
theory for Heisenberg spins. This procedure leads to a Tc (lambda) /Tc(0) almost 
identical to that shown in Fig. 2; however the projection method Tc is lower than 
those shown in Fig. 2. For example, the direct-integration Tc at A = /i = is 
O.lt, much less than the 0.17t shown in Fig. 2. A numerical error originally led us to 
believe the Tc's of the two approaches coincided. The discrepancy may most easily 
be understood by expanding F = — Tln[ZLoc] to order al^. The result is a quadratic 




(25) 




12 



form(5F ~ EmnainAnmaim- Foi example, at g = 0, Amn = (5mn(l-l/3(aon)^) + 3^^^^. 
Tc is the temperature at which A first has a zero eigenvalue. The prjection method 
result corresponds to setting = '^mn(l — l/3(aon)^) and ain = constant; in other 
words it produces a lower Tc because it does not permit an optimal choice of ain- 
We have therefore not used the direct integration method in this paper. We note, 
however, that the projection method provides a transparent and physically appealing 
motivation for the result, found also in the detailed calculations presented below, that 
the Tc is controlled by the kinetic energy at Tc, is borne out by our detailed solution 
of the model. 



In this section we discuss the qualitative behavior of the solutions of Eqs. |T6[ 
Much of the behavior is similar to that found in I. The new feature is the physics 
of double exchange, which is expressed via bi, via the angular integral and via the 



At T ^ 0, the 9 integral is dominated by the regime cos^ = 1, so bi = bg — (u^+yu). 
From Eq. |19] one sees that at T = the antialigned component of G vanishes, while 
the aligned component is determined by (bo + bi) which is given by an equation 
identical to that considered in I. Therefore, all of the results obtained in I for the 
T ^ limit hold also here. At T > Tc, there is no long range magnetic order. Thus 
bi = 0, there is no ^-dependence and bo is given by an equation which differs by a 
factor of I from that treated in I. 

Further insight into the quantity bi may be gained from the A = limit. At T = 
and A = the quantity bo + bi is found from Eqs. [1^ to be 



This is precisely the usual non-interacting solution: ImGLoc 7^ in a semicircular 



III. QUALITATIVE BEHAVIOR 



factors of \ on the right hand side of Eqs. |TB[ 




(26) 
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band of full width 4t. In the present conventions, the fermi level is at = and for 
A = the maximum of the spectral function is at a; = — /x. The self-energy for this 
solution vanishes. 

At T > Tc and A = 0, bi = and 



+ /i - iJ2 - (u; + ii)' 



(27) 



Here ImGLoc 7^ in a semicircular band of full width 1\p2X: the fact that neighboring 
spins are uncorrelated has reduced the bandwidth, and thus the kinetic energy, by a 
factor of \f2. This may also be seen by a direct evaluation of K from Eq. |23|. Further, 
the self energy is 



S(^) = -bo = - ^2 - (^ + - -(^ + /x) (28) 

and has a non-zero imaginary part at the fermi surface (cu = 0), corresponding physi- 
cally to scattering by spin disorder. However, this scattering is not too strong. From 
Eqs. |2^, [TT] one finds that the product of the imaginary part of the self energy and 
the density of states at the Fermi level is (2 — /i^)/vr. This number is rather less 1, 
and implies a mean free path longer than pp^. This spin disorder scattering decreases 
as T is decreased below Tc. 

The model with A = was studied in the dynamical mean field method by Fu- 



rukawa [jT5[, who obtained Eq. |2^. Furukawa also used a method he referred to as 
solving the equations at constant magnetization to produce an interpolation formula 
describing the temperature dependence of E" for < T < Tc. We believe these re- 
sults are similar but not quite equivalent to those we obtain by solving the equations 



T6| directly. However, the minor differences between Furukawa's results and ours are 
not important. The main point is that the scattering at T > Tc predicted by this 
calculation is much too small to explain the data. 

One may calculate Tc at zero coupling by linearizing the second of Eqs. in bi. 
One finds that Tc(/i) is given by the solution of 
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where f is the fermi function. 

We now return to the issue of the effects of the electron-phonon couphng. At T = 
the mean field equation is identical to that considered in I. Prom this work we learn 
that there arc three regimes: weak coupling, in which limT^o^^(T) = 0, limT->oP(T) = 
and dp/dT|T=o ~ A; intermediate coupling, where < limT_»ox^(T) < ~ 1, 
< limT_»op(T) < oo and dp/dT|T=o may have either sign, and strong coupling, 
where < limx^ox^(T) and limT^oP(T) = oo. Here Xc is the value of frozen in 
lattice distortion above which a gap appears in the electron spectral function. In the 
strong coupling regime one may think of the electrons as being localized as polarons. 

Another crucial result of I is that the transitions between the different regimes 
are controlled by the values of an effective coupling determined by the ratio of an 
electron-phonon energy to a kinetic energy. As we have seen, the kinetic energy is 
temperature dependent because of double exchange; thus as temperature is varied 
the behavior of the model may change from "metallic" (dp/dT > 0) to insulating 
(dp/dT < 0). As T is decreased below Tc there are two effects causing a decrease 
in the resistivity: the spin scattering freezes out and the effective electron-phonon 
coupling weakens. 

IV. HALF-FILLING 

In the section we present and discuss results of numerical calculations for the 
particle-hole symmetric (n = 1) case. We begin with ferromagnetic Tc shown in Fig. 
2. One sees that Tc decreases with increasing A; the variation is particularly rapid in 
the region A ~ 1 which is shown below to be the critical value at which the model 
goes from metal to insulator. 

For n = 1 and all A we verified that the transitions were second order by compar- 
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ing the Tc obtained in this manner to the Tc obtained by determing the temperature 
at which the non-magnetic state becomes hnearly unstable. We also checked for 
metastability at various n and A by starting our iterations with saturated magneti- 
zation (bi = bo — (cu + /i)) and with very small bi, and verifying that both initial 
conditions converged to the same solution. The magnetic transition was always found 
to be second order. 

It is interesting to compare our results to those of Ref. [^]. The method used 



by these authors to treat the magnetic fluctuations is very similar to the "projection 
method" discussed at the end of section II. We found that the method did not give an 



accurate value for Tc but did reproduce the coupling dependence well. Ref p3| used a 
model with one orbital per site; we should therefore compare their results for n = 1/2 
to ours for n = 1. Their quantity Cp = Ajt/2K corresponds precisely to our A; the 
factor of two comes from the orbital degeneracy as explained in our previous paper 
I. As far as can be determined from Fig 1 of ref ||23|, their calculated Tc(ep)/Tc(0) 
agrees very well with our Tc(A) /Tc(0). The correspondence is interesting because the 
calculation of ref [^] was done with quantum phonons with the rather high frequency 
oj = 0.5 (units not specified, but presumably set by the electron hopping, t). This 
supports our claim that quantum effects are not important at temperatures of the 
order of Tc. We believe the results presented in p3[ for the doping dependence of Tc 
are not physically relevant because they are based on a model with one orbital per 
site, which therefore has no kinetic energy at n = 1. They argue that the one-orbital 
model is justified by the existence of the Jahn- Teller splitting. Our results show that 
this is not the case. We discuss the physics of the doping dependence of Tc further 
in the conclusion. 

We now turn to the temperature dependence of the resistivity, shown in Fig. 
3. The curves display kinks at the ferromagnetic Tc. The resistivity drops as T 
is decreased below Tc both because the magnetic contribution to the scattering be- 
gins to decrease at Tc and because the effective electron-phonon interaction becomes 
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weaker. Prom these curves we may distinguish "metaUic" (dp/dT > 0) and "insulat- 
ing" (dp/dT < 0) regimes. At T > Tc A = 1 marks the boundary between metaUic 
and insulating regimes; at A = 1, dp/dT = 0. For T < Tc the crossover occurs at 
the somewhat larger A ~ 1.15. The difference in the critical A required to produce 
insulating behavior reflects the effect of spin alignment on the electron kinetic energy. 
We also note that although it is difficult to perceive on the logarithmic scale used 
in Fig. 3, for 1.08 < A < 1.15, limT-+op(T) = Po is neither zero nor infinite. For A 
sufficiently close to 1.15, hmx-^o^ < 0. Similar behavior was discussed at length in 
I. 

Figure 4 shows the temperature dependence of the mean square lattice distortion 
x^. From this one may distinguish the low-T weak, intermediate and strong coupling 
regimes, based on the T — > limit of r(T). The regimes were discussed at length in 
I. Roughly, in weak coupling r(T = 0) = 0; in intermediate coupling < r(T = 0) < 1 
and in strong coupling r(T = 0) > 1. In intermediate coupling there is a frozen-in 
lattice distortion which affects the T = physics but is not large enough to open a gap; 
for strong coupling the distortion is large enough to open a gap and cause insulating 
behavior. For r^ > 0.25, hmx-^odp/dT < even though if r^ < 1, limT^oP(T) is 
finite. 

The effects of double exchange may be seen in Fig. 4. Tc is visible as a kink 
on each curve. At T > Tc, dr^/dT decreases, implying a stronger electron-phonon 
coupling. For A > 0.9 the T = values obtained by extrapolating the T > Tc 
curves to are non-zero, and are higher than the actual T = values, because the 
reduction of kinetic energy due to spin disorder has effectively made the electron- 
phonon coupling stronger. Note that the T = extrapolation of T > Tc portion of 
the curve corresponding to A = 1.05 is about r^ = 0.95. This is slightly less than 
the critical value — i found in I to mark the boundary between finite and infinite 
p(T = 0) at n = 1. From this we would infer at T > Tc, dp/dT changes sign at 
A > 1.05, as indeed is seen in Fig. 3. 
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The curves presented in Fig. 4 show r^ in arbitrary units. To estimate the magni- 
tude of the effect in Rei-xAxMnOs we note that in LaMnOs each 0-ion is displaced 
^ 0.15A from its ideal perovskite position [|^. The estimates obtained in |^ imply 
A ~ 1.3 — 1.5 in that material; thus r^ = 3 in fig. 4 corresponds to an rms displacement 
of an 0-ion of about O.ISA. 

We next consider the temperature dependence of the kinetic energy shown in Fig. 
5. At A = the kinetic energy changes by about l/\/2 = 30% between T = and 
T = Tc, and has a weak T dependence at T > Tc. For A = 0.71 the kinetic energy 
changes between and Tc by a somewhat larger amount; for A = 1.11, by still larger 
amount, for A = 1.29, yet larger. These changes come from the previously discussed 
interplay between double exchange and electron-phonon coupling. As T is increased 
from zero, the spins disorder. This reduces the electron kinetic energy and permits 
the electron-phonon coupling to further localize the electrons, reducing their coupling 
yet more, etc. We also note that we found the ratio between T^. and the kinetic energy 
at Tc to be the same within a few percent for all n and A studied. For n = 1 this can 
be seen by comparing Figs 2 and 5. 

These arguments also explain the magnetic field dependence of the resistivity. 
Increasing the field aligns the spins, increases the kinetic energy, and decreases the 
effective electron-phonon coupling, leading to a large change in resistance as shown 
in Figs 6a and 6b. When the decrease in effective electron-phonon coupling tunes the 
model across the "metal" -"insulator" transition, as in Fig. 6b, the magnetoresistance 
is particularly large. 

Further insight into the interplay of double exchange and localization comes from 
the optical conductivity shown in Fig. 7. Panel 7a shows <y{oj) at different T for 
the weak coupling A = 0.7. At low T a has approximately the Drude form cr(co') = 
r/(u;^ + r^) with scattering rate F ~ T as expected from classical phonons. As 
T is increased through Tc the Drude peak broadens and acquires a T-independent 
part, due to spin scattering. Panel 7b shows a for the moderate coupling A = 1. 
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At low T, cr has the Drude form; as T is increased a broad peak develops; this is 
due to transitions between the two Jahn- Teller split levels. It is broad because the 
phonon coordinate is strongly fluctuating, so the level position is not well deflned. As 
T increases beyond Tc the peak broadens almost to indistinguishability. Note also 
that as T is increased, the optical spectral weight / dua^uj) decreases, reflecting the 
increasing localization of electrons by phonons. In models such as the present one 
which do not have Galilean invariance and involve only a limited number of orbitals, 
the f-sum rule spectral weight is not constant and is indeed proportional to the kinetic 



energy ||2^. Panel 7c shows that at a stronger coupling a does not have the Drude 
form, and the peak is already evident at Tc/2. Note that the maximum in a has 
moved to a slightly higher frequency. Recently published data of Okimoto et al. on 
Lai 825Sr lysMnOs (H are similar to the curves shown in Fig. 7c, although our use of 
classical phonons means that we cannot obtain the very narrow Drude peak found at 
low T. Panel 7d shows a at the still stronger coupling A = 1.15 where the model has a 
large frozen-in lattice distortion even at T = 0. The a has an insulating appearance, 
above and below T^,, but as T is decreased the peak in a shifts to a lower freqency 
and grows in intensity, reflecting the effectively weaker coupling. The nonmonotonic 
behavior of the dc conductivity is not reflected in (j{uj) at uj > 0.5. The curves in Fig 



7d resemble data recently obtained by Kaplan et. al. on Ndo.TCao.sMnOa . 

We will discuss the physical interpretation of cr{uj, T) in more detail in the next 
section and in the conclusion. Here we note that in the strong-coupling regime the two 
d-states on a site are split. The peak in the optical conductivity corresponds roughly 
to a transition in which an electron moves from an occupied orbital on one site to 
an unoccupied orbital on an adjacent site. In our classical approximation, "Franck- 
Condon" transitions involving also emission or absorption of a phonon cannot occur at 
all. In a more realistic model such effects would e.g. increase the low-frequency tails 
by a small amount. The width of the peak in (j{uj) is determined by the broadening of 
the localized states due to electron hopping and by thermal broadening, which leads 
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to a range of lattice distortions and thus to a range of splittings. 

V. DIFFERENT DOPINGS 

In this section we present and discuss results of numerical calculations for the 
particle-hole asymmetric case n 7^ 1. As discussed in I (see especially Figs. 10, 
11), at n 7^ 1 in the strong coupling limit the spectral function has a three-peaked 
structure. The outer two peaks represent the Jahn- Teller-split Cg levels on occupied 
sites, and occur also for n = 1. The middle peak comes from unoccupied sites, on 
which there is no Jahn- Teller splitting. These states tend to fill in the gap created 
by the Jahn- Teller splitting and mean that stronger coupling is required to obtain 
insulating behavior at n 7^ 1 than at n = 1. Further, in the strong coupling limit the 
temperature dependence of physical quantities is determined by the energy difference 
between filled and mid-gap states; thus at fixed Jahn-Teller splitting the activation 
gap for physical properties is much less at n 7^ 1 than at n = 1. Note also that T^ 
is controlled by the electron kinetic energy which is in turn controlled by the Jahn- 
Teller sphtting. Therefore, in the strong coupling limit at fixed Tc the activation gap 
characterizing the T > Tc resistivity is much larger at n = 1 than at n 7^ 1. 

This physics is immediately apparent in the resistivity curves for n = 0.75 and 
n = 0.5 shown in Fig. 8. Comparison to Fig. 3 shows that much stronger couplings 
are required to obtain a dp/dT < for n = 0.75 than for n = 1 and stronger couphngs 
yet are required for n = 0.5. The smaller value of the activation gap relative to Tc 
means that the resistivity rises less before the behavior changes at Tc for n 7^ 1 than 
for n — 1. For n = l, A = l.lwe found an order-of- magnitude rise in p as T is 
decreased to Tc, and we found metallic behavior below Tc. At n 7^ 1 it is difficult 
to produce much of an up-turn in p at T > Tc for parameters such that the model 
is metallic at T = 0. This physics is due to the particular (Jahn-Teller) form of the 
electron-phonon coupling we have chosen to study. 
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The same arguments mean that it is not possible to get as large a magnetore- 
sistance at n 7^ 1 as at n = 1. Figure 9 shows the temperature dependence of the 
resistance for several different coupling strengths and magnetic fields. To get even a 
moderately large effect one must choose a very strong coupling, such that the model 
is insulating for both T > Tc and T < Tc- 

Finally, Fig. 10 displays the temperature dependence of the optical conductivity 
at n = 0.75 and moderate (A = 1.29) and strong (A = 1.49) coupling, and compares 
this to the momentum-integrated spectral functions. One sees by comparing ener- 
gies that at strong coupling the two maxima in the conductivity may be associated 
with transitions from the lowest peak in the spectral function (representing occupied 
orbitals on occupied sites) to the middle feature (representing unoccupied orbitals 
on unoccupied sites) and to the higher feature (representing unoccupied orbitals on 
occupied sites). At n = 1 the middle peak in A is absent and a has only one peak, 
as seen in Fig. 7d. Of course, the on-site d-d transition is not optically active: the 
calculated conductivity involves electron motion from one site to another, and for this 
reason a is not simply given by a convolution of two local spectral functions. One 
may see this in Fig. 10. The central peak in kiio) has less area than the upper one, 
yet the lower peak in the corresponding cr(a;) is the larger. This may be understood 
from the above arguments: a transition from the lower to the middle peak of A (a;) 
necessarily involves moving an electron from one site to the other, but some of the 
transitions from lower to higher are on-site transitions which do not contribute to a. 

Note also that the T = 2Tc spectral function shown in Figs 10b,d has a sharp 
minimum at a; -|- // = 0. This is a consequence of the fact, discussed in I, that the 
probability of a small-amplitude lattice distortion is small because of the xdx measure, 
and decreases as T increases, due to the shift to higher (x^) of P(x). As can be seen 
in the corresponding optical conductivity curves in Figs 10a,c, this minimum is of 
little significance for other physical quantities. 
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VI. CONCLUSION 



We have used the "dynamical mean field" approximation to solve a model of elec- 
trons ferromagnetically coupled to classical spins and Jahn- Teller coupled to localized 
classical oscillators. In a companion paper (I) we considered electron-phonon coupling 
in a variety of models without double exchange. The results presented in Section IV 
for the half-filled case bear a striking resemblance to data for the "colossal magnetore- 
sistance" materials Rei-xAxMnOs in the 0.2 < x < 0.5 regime where the ground state 
is metallic. We believe the agreement supports the idea that the important physics 
of Rei-xAxMnOs involves the interplay between a strong electron-phonon coupling 
and the "double exchange" effect of magnetic order on the electronic kinetic energy. 
Specifically, the p(T) curves shown in Fig. 3 are very similar to those shown, e.g in 
||29| , |30[| . Varying the electron-phonon coupling produces changes very similar to those 
found experimentally by varying x and the constituents Re and A. The magnetic field 
dependences shown in Fig 6 also bear a striking resemblance to data. Fig 6b looks 
very like Fig. 2 of [Q, while Fig 6a resembles magnetoresistance data which would 
be observed for Lai.6Sro.4Mn03. (It should be noted, however, that the fields used to 
produce our curve, although very small compared to microscopic energies, are larger 
than experimental fields by a factor of about 5). The variation of the rms lattice 
distortion shown in Fig 4 has been observed via measurements of the Cg component 
of the oxygen Debye- Waller factor []3ll , ^ |. 

Further, optical conductivity data of Okimoto et. al. |^ on Lai.825Sr.i75Mn03 
bear a strong qualitative resemblance to Fig 7b, while data obtained by Kaplan et. 



al. p8[ strongly resemble Fig 7c. As noted in sections IV and V, in our interpretation 
the higher frequency peak in (j{uj) is due to transitions between levels split by an 
electron-phonon coupling. Okimoto et. al. interpreted the higher peak differently, 
attributing it to transitions from an inital Cg state aligned to the core spin to a 
final eg state antialigned to the core spin. They argued that their identification was 
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supported by the fact that in their data the higher peak was only visible at T > Tc, 
and vanished at low T when all spins were aligned. However, the data of Kaplan 
et. al. demonstrate that in some samples the high frequency peak does not vanish 
below Tc and indeed grows in oscillator strength as T decreases. This rules out the 
interpretation of Okimoto et. al., at least for that sample. 

The detailed qualitative agreement between data and our model leaves little doubt 
that we have identified the important physics governing the Rei-xAxMnOa materials. 
However, several very important issues remain unresolved. One concerns the origin of 
the experimentally observed material and doping dependence of the results, which are 
modelled in the n = 1 calculations by varying the electron-phonon coupling. Another 
is the degree of "fine-tuning" of parameters required. A third concerns the effects 
of omitted interactions and a fourth is that, as shown in Section V, computations 
at different electron concentrations n 7^ 1 agree much less well with data. In the 
remainder of this paper we present a qualitative discussion of all of these issues, 
which, we argue, are closely related. 

We begin with the n 7^ 1 calculations. We showed that the differences between the 
n = 1 and n = 0.75, 0.5 results are due to the presence, for n 7^ 1, of mid-gap states 
in the spectral function (shown, e.g., in Fig. 10). These mid-gap states occur because 
we used a particular form of electron-phonon coupling, namely a Jahn- Teller coupling 
which splits the d-state degeneracy on a site if there is one electron on the site, and 
does nothing otherwise. In ref [^] it is argued that the existence of the Jahn- Teller 



coupling justifies a model involving only one orbital per site. The results presented 
here suggest that this is oversimplified, because it does not take in to account the 
mid-gap states. A model with only Jahn- Teller coupling does not suffice. However, 
results presented in I strongly suggest that if the model were extended in a way 
which moved both the upper peak and the mid-gap states up in energy, then the 
model would become effectively a single-orbital model and results for n 7^ 1 would 
much more closely resemble those obtained for n = 1. 
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One omitted piece of physics which will have precisely this effect is the breathing- 
mode distortion of the oxygen octahedron around an Mn site. The breathing mode 
couples to charge fluctuations on the Mn site. This coupling is likely to be at least 
as strong as the Jahn- Teller coupling, as may be seen from the following argument: 
the Jahn- Teller coupling is due to the dependence of the force exerted on an O ion on 
the orbital occupied by the outer-shell d-electron. Whatever its magnitude, this force 
is unlikely to be larger than the force created by simply removing that d-electron, 
and making an unbalanced charge. The breathing-mode coupling was recently argued 
to be important for the small-x structural phase boundary |20|. To understand the 
effects of the breathing mode, consider again Fig. 10. The central peak in the spectral 
function depicted in the low-T curve in Fig lOd gives the states available for adding 
an electron to an unoccupied site. If such a site has a breathing distortion already 
present, the energy cost of adding an electron will be increased, thus the middle 
feature will also move up in energy, increasing the gap as required. 

Another important piece of physics is the on-site Coulomb interaction. This must 
be strong because if it were not, the Hund's coupling Jh would not be large [|T7|. 
The Coulomb interaction leads to two related effects. One is most easily discussed by 
reference to the spectral functions and optical conductivities shown in Fig. 10. Now 
the a; > part of the spectral function corresponds to states into which an electron 
maybe added; the upper peak thus gives the states available for adding an electron 
onto a site which already has an electron. The Coulomb interaction must move such 
states up in energy, and must similarly move up the second peak in a{uj). If the 
Coulomb energy is of the order of the Hunds coupling, then it is very likely that this 
effect will move the higher peak out beyond the physically interesting energy range 
uj <3 eV. 

The combined effect of the breathing distortion and the Coulomb interaction is 
therefore to lead to a spectral function with at most two peaks in the energy range 
of interest. The only difference between this realistic situation and the situation 
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encountered in the n = 1 calculations is that the realistic spectral function is not 
symmetric under the interchange of the two peaks. This asymmetry was shown in I 
not to be important. 

A second effect of a strong Coulomb interaction is to localize the electrons. It is 
likely that the observed very strongly insulating behavior of ReMnOs is not due solely 
to the Jahn- Teller order, and that ReMnOs is to some degree a Mott insulator. Now 



the kinetic energy K of a Mott insulator has a pronounced doping dependence |£7 
For Rei_xAxMn03 one would expect K(x) to increase with x for x < 0.5. Because, as 
we have argued at length, the properties of electron-phonon models are controlled by 
the ratio of a coupling energy and a kinetic energy, this will lead to an x-dependence 
of the effective coupling strength, with larger x having a weaker effective coupling. 
We believe that this strong x-dependence of the effective coupling accounts for the 
ubiquity of the "colossal" magnetoresistance phenomenon. Different materials have 
different bare electron hoppings and probably different electron-phonon couplings, 
but in all materials the variation of the electron kinetic energy with x is large enough 
to sweep the effective coupling through the critical value at some x between 0.1 and 
0.5. 

The breathing distortion may be studied via the dynamical mean field formalism 
used here; one must simply integrate over another variable in Eq |T2|. The on-site 
Coulomb interaction may also be included, but one must perform functional integrals 
rather than simple integrals. Monte Carlo techniques are required, the computational 
expense is greater and the accuracy is less. Such an investigation would however be 
desirable. 

Two other effects not included in the calculation should be mentioned. Quantum 
fluctuations of the phonons have been omitted. As discussed in I, these will in the 
absence of long range order or commensurate density lead to metallic behavior at 
sufficiently low T, even in the strong coupling limit. The neglect of the phonon 
momentum and quantum fluctuations of core spins is not an important approximation 
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because we are primarily interested in phenomena at temperatures of order room 
temperature; however if needed they could be incorporated into the formalism. The 
neglect of intersite phonon correlations is potentially more serious. It is tempting to 
argue that they are unimportant because we are interested in optical phonons, which 
are usually weakly dispersing. However, in the ReMnOa structure each O is shared 
by two Mn; there must thus be a strong correlation between Jahn- Teller distortions 
on adjacent sites. In LaMnOs the Jahn- Teller distortions have long range order; 



estimates presented in pO| suggest that in Rei-xAxMnOs the correlation length of the 
Jahn- Teller distortions is ~ as long as the resistivity is well above the Mott limit. 
Extending the present calculations to include the effects of intersite correlations is an 
important open problem. It is worth addressing because in the present calculations 
the correlation length is zero and the strong-coupling physics is of polarons. In the 
infinite correlation length limit, the physics has to do with interband transitions in a 
bandstructure defined by Jahn- Teller order. The situation in the actual materials is 
presumably intermediate between these two limits. 
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APPENDIX A: ANALYSIS OF OBSERVED RESISTIVITY 



In this Appendix we discuss the observed T > Tc resistivities of Rei-xAxMnOs. 
We note that the observed strong x-dependence suggests that the number of active 
carriers is x. For x classical particles hopping with probability W on a cubic lattice of 
lattice constant a, 

^ = 3^ 

Using a ~ 4A as appropriate to Rei-xAxMnOs we have 

Kb i xp{\l — cm) 

/^From this equation one may easily see that observed resistivities, which are typ- 
ically greater than 0.01 Q — cm at T > Tc and x < 0.3, and increase rapidly with 
decreasing x imply values of hW /k-^T much less then unity. If ^iW/keT ^ 1, a 
particle has time to thermalize before it moves, and a classical model is appropriate. 



APPENDIX B: RESISTIVITY OF DOUBLE EXCHANGE ONLY MODEL 

In this Appendix we consider in more detail the resistivity of the double exchange 
only model. In this model, resistivity comes from spin disorder. It is maximal at 
T ^ Tc and vanishes at T = 0. The resistivity has been calculated, using methods 
which are essentially perturbative in the amplitude of the spin disorder, by Kubo and 



Ohata |TJ] and more recently by two of us and Littlewood [T^. The spin scattering 
was found not to be too strong. As discussed in the text, similar results have been 
obtained using the dynamical mean field method. These calculations omit several 



physical effects and have been questioned recently by Varma |T^. In this Appendix 
we show that the omitted effects are not important. 

We begin by describing the omitted effects. In the double exchange only model, the 
scattering is due to spin disorder which, if the core spins are assumed to be classical, 
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may be treated as static scattering with the important proviso that the disorder 
is annealed, not quenched. When apphed to a model with static scattering, the 
dynamical mean field approximation with semicircular density of states is equivalent 
to the coherent potential approximation (CPA) for the Bethe lattice |21|. The CPA 
neglects localization (as do the perturbative calculations [0,^). The lack of closed 



loops on the Bethe lattice also means that Berry phase effects arising from particle 
motion in a spin background are omitted. 

We consider the Berry phase effects first. In the double exchange model the 
hopping matrix element between two sites i and j is with core spins characterized by 
polar angles (^i, 0i), (6'i, 0j), is 

tij = t(cos|cos^ + sin^sin^e'(*-'^^)) (Bl) 

If closed loops are not important one may choose the ^\ independently on each site 
and recover the familiar double exchange result ty = tcos(^^i — Q-^jl. In general 
the factors around a closed loop produce something like a magnetic field, which 
may scatter electrons. In the limit of strong ferromagnetic correlations the phase 
dependent term may be seen to be very small because all nearby sites have very 
similar angles, which may be taken to be near 0. In the limit of uncorrelated spins we 
may estimate the size of the effective field by comparing the phase sensitive part of the 
hopping to the phase insensitive cos(0i — 0j)/2 part. By integrating ty, around square 
placquette one finds that the phase sensitive part is |^e^"^' while the phase insensitive 
part is t^/4. Thus the rms deviation of the amplitude for an electron to move around 
a placquette is l/4-\/2 of the phase insensitive part. This combined with the relative 
insensitivity of three dimensional physics to closed loops suggest phase effects, while 
interesting, are too weak to cause the observed strongly insulating behavior. 

We now turn to localization. The problem at hand concerns electrons with random 



hopping, which has not received much attention. Economou and Antoniou p3 have 



studied a Bethe-lattice model in which the hopping amplitude t has the symmetrical 
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distribution 



2 



PE(t) = — Vt?-(t-to)2 (B2) 

For this model t, the mean value of t, equals to and the variance ((t — t)^) = tf/4. 
The double exchange model at T >> Tc (so the spins are completely disordered) 
corresponds to the distribution 

2t 

Pd-cx(t) = ;^^^(tD - t) (B3) 

The localization effects of the double exchange distribution have not been determined. 
We expect that because the most probable value is also the largest hopping, the 
double exchange distribution will produce fewer localized states than a semicircular 
distribution with the same mean and variance. Now from Eq ^3| one sees that the 
double exchange distribution has mean t = 2tD/3 and variance of t^/18. Thus it 
should produce fewer localized states than the model of Economou and Antoniou 
with to = 2tD/3 and ti = v^tn/B, i.e. with ti/to = 1/v^. Inspection of Ref pi | 



reveals that at this ratio of ti/to, a negligible fraction of the states are localized. 
We therefore conclude that localization effects are not important. Ref on the 



contrary asserts that the double exchange model with completely disordered spins is 
better modelled by the Economou- Antoniou distribution with ti/to somewhat larger 
than unity, so a non-negligible fraction of the states are localized. Drawing precise 
conclusions is somewhat difficult because one result of Ref is that the number of 
localized states increases rapidly for ti/to > 1. Nevertheless, we believe the estimate 
ti/to ~ 0.7 obtained above shows that localization effects are unlikley to be important. 
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APPENDIX: FIGURE CAPTIONS 



Figure 1. Qualitative temperature (T) - doping (x) phase digram of Rei_xAxMn03, 
with magnetic phases (F =ferromagnet, AF=antiferromagnet, P=paramagnet), 
structural phases (JT =Jahn- Teller order, no label=no order), and transport 
regimes (M = "metal", dp/dT > 0, I = "insulator" , dp/dT < 0) indicated. 
The sohd hues are magnetic phase boundaries, the heavy dashed hne is the 
Jahn- Teller boundary and the light dotted line is the metal insulator crossover. 
For X > 0.5 different physics, involving charge ordering, is important at low T. 
Different materials may have phase diagrams differing in some details, and the 
magnetic and structural boundaries may not coincide at low T. 

Figure 2. Dependence of ferromagnetic Tc on coupling constant for n = 1 (heavy solid 
line), n = 0.75 (hght sohd hne), n = 0.5 (light dashed line). The analytic zero- 
coupling results are indicated by dots; the analytic strong coupling Tc = 
results by the heavy dotted line for n = 1. Only for n = 1 do the numerical 
calculations extend into the strong coupling regime. 

Figure 3. Temperature dependence of resistivity at n = 1 for couplings A = 0.32 (lowest 
curve), 0.71, 1, 1.08, 1.12, 1.15, 1.20 (highest curve). 

Figure 4. Temperature-dependence of mean-square lattice distortion for n = 1 and cou- 
plings A = 0.71 (lowest), 0.9, 1.05, 1.12, 1.2 (highest). 

Figure 5. Temperature (T) dependence of electron kinetic energy (K) for n = 1 and 
A = (second lowest curve), 0.71, 1.12 and 1.29 (highest curve). The lowest 
curve corresponds to A = in the model without double exchange. 

Figure 6. Temperature dependence of resistivity at different values of magnetic field, h, 
for A = 0.7 (Figure 6a) and A = 1.12 (Figure 6b). The parameter h is related 
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to the physical field hphys by h = gA*BSchphys/t. Using g = 2, t = .6 eV and 
Sc = 3/2 means h = 0.01 corresponds to hphys = 15T. 

Figure 7. Optical conductivity, n = 1, T = 0.02 (light solid line), T = Tc/2 (light dashed 
line), T = 3Tc/4 (light dotted line), T = Tc (heavy solid line), Tc = 2Tc (heavy 
dashed line). Panel a: A = 0.71 (T^ = .15), panel b: A = 1 (T^ = 0.10), panel 
c: A = 1.08 (Te = 0.08), panel d: A = 1.15 (Te = 0.0675). Note that in panel a 
the lowest T is .025 not .02, and cr(a; = 0) for this curve is 21.4. 

Figure 8. Resistivity (p) vs temperture (T) for n = .75 (upper panel) and n = 0.5 (lower 
panel) and couplings A = 0.71 (lowest), 1.12, 1.41, 1.49, 1.58 (highest). 

Figure 9. Magnetic field dependence of resistivity for n = 0.75 and A = 1.12 (lower panel), 
A = 1.46 (middle panel) and A = 1.49 (upper panel). 

Figure 10. Optical conductivities and spectral functions for n = 0.75, A — 1.29 (panels 
a,b) and A = 1.49 (panels c,d). Panel a: a{u;), A = 1.29 and T = .04 (fight 
solid line), T = .061 (light dashed line), T = .081 = (light dotted line) and 
T = .162 (heavy solid line). Panel b: spectral function, A = 1.29 and T = 0.04 
(solid fine) and T = 0.162 (dashed fine). Panel c: a{uj), A = 1.49 and T = 0.02 
(light solid line), T = .028 (light dashed line) T = .045 (light dotted line), and 
T — 0.059 = Tc (heavy solid line) and T = .115 (heavy dashed line). Panel d: 
spectral function, A = 1.49, T = 0.02 (solid fine) and T = .115 (dashed line). 
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